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Abstract 

We propose a simple mechanism which enforces \Uaj\ = \U T j\ Vj = 1, 2, 3 in the 
lepton mixing matrix U. This implies maximal atmospheric neutrino mixing and a 
maximal CP-violating phase but does not constrain the reactor mixing angle #13. 
We implement the proposed mechanism in two renormalizable seesaw models which 
have features strongly resembling those of models based on a flavour symmetry group 
A(27). Among the predictions of the models, there is a determination, although 
ambiguous, of the absolute neutrino mass scale, and a stringent correlation between 
the absolute neutrino mass scale and the effective Majorana mass in neutrinoless 
double-beta decay. 
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1 Introduction 



With the recent results of the Double Chooz, Daya Bay, and RENO Collaborations pQ 
the earlier hints [2] of a non-zero reactor mixing angle #13 have been confirmed. The 
unexpectedly large value of #13 pfl EJ 0] renders a /x-r interchange symmetry in the neutrino 
mass matrix [5], and therefore also tri-bimaximal mixing j6], highly unlikely, since that 
symmetry was tailored to achieve #13 = at some energy scale. However, there is a 
different version of the /i-r interchange symmetry, which is based on a generalized CP 
transformation that includes the \l—t interchange [3, [8]. In this version, the maximal 
atmospheric mixing angle #23 is not coupled with a vanishing #13 but rather with a maximal 
CP-violating phase S in the lepton mixing matrix. Phenomenologically, this scenario is 
fully viable. 

In this letter we introduce a new mechanism for generating this type of lepton mixing. 
In order to establish our notation we firstly define the lepton mass Lagrangian as 

Anass = -4MA + ~ ulC~ X M v U L + H.C., (1) 

with Mi and M. v being the mass matrices of the charged leptons and of the neutrinos, 
respectively; the latter mass matrix is of the Majorana type. Those mass matrices are 
diagonalized by 3 x 3 unitary matrices Ut and U v according to 

U}M e M}U e = diag (m 2 e ,m J, m 2 T ) , (2a) 

UjMMu = diag (mi, m 2 ,m 3 ) , (2b) 
respectively. Then the lepton mixing matrix U is given by 

U = U\U V . (3) 

Our idea is the following. Suppose that we have a model which gives a real matrix U v 
and 

i I 1 1 M 

U t = = -= 1 u u 2 }, (4) 
V3 \ 1 co 2 u J 

where to = (— 1 + i\/3) / 2. Then, it is trivial to see that the mixing matrix U has the 
property 

\Uri\ = \U rj \, Vj = 1,2,3. (5) 

When using the standard parameterization of the mixing matrix [9], the relations (151) 
require [3, [8] 

cos #23 = sin 9 2 3 = -^=, (6a) 

v 2 

sin #i3 cos 5 = 0, (6b) 

whence it follows 

O23 = \, (7a) 
5 = ±|, (7b) 
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since we know that #13 ^ 0. 

This paper is organized as follows. In section [2] we develop two seesaw models based 
on the idea laid out above. These models make predictions beyond those in (J7|); the extra 
predictions are presented in section [3j The conclusions of the paper are summarized in 
section HI 



2 The models 

The fermion sectors of our models contain the usual leptonic Standard Model multiplets, 
namely three left-handed gauge- SU (2) doublets, subsumed under the symbol Dl, and 
three right-handed charged-lepton gauge singlets, subsumed under the symbol £ R . The 
scalar sectors contain three Higgs doublets with weak hypercharge 1/2, which we subsume 
under the symbol 0. For the symmetry transformations of the models we make use of the 
matrices 
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We define the bilinears 

Wr] q = Mm + $4ir + fotsR, (9a) 

= faiw + wfatui + u 2 fohR, (9b) 

Wr] 2 = <Pih R + ui 2 <f) 2 e 2R + wfaltR, (9c) 

which are analogous to the ones used in models based on the symmetry A4 [10]. Indeed, 
under the transformations 

S : £ R ->• A£ R , -»• A<f>, (10a) 

T : l R -> E£ R , ->■ E(f), (10b) 



'11a) 



the bilinears transform as 

Wr], With, 

[0^]. A ^[0^] r (lib) 

These transformation properties allow us to write down the charged-lepton Yukawa La- 
grangian 

3 

3=1 

if we supplement (110b[) by 

T : Dl^&Dl. (13) 

The charged-lepton Yukawa Lagrangian (fT2l looks very similar to the one of some A4- 
based models [El ITT] , but that look is misleading — in our models the roles of D L and £ R 
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are reversed relative to the A 4 -based models. This can also be seen by computing the 
mass matrix of the charged leptons, which in our models is 

M e = diag (hi, h 2 , h 3 ) (y3uSj diag (v x ,v 2 , v 3 ) , (14) 

where Vj denotes the vacuum expectation value (VEV) of the neutral component of the 
Higgs doublet <pj. Thus, U^M^ is diagonal if hi = h 2 = h 3 , i.e. if the Yukawa coupling 
constants are all equal, whereas in the A A models one needs equality of VEVs. The 
equality of the hj is achieved by assuming invariance of C Y under 

T' : D L ^ ED L , £ R -> Ci R . (15) 

Therefore, in the following we shall use 

hi = h 2 = h 3 = h, (16a) 
TJlMt = V3hdia,g(vi,v 2 ,v 3 ). (16b) 

We emphasize that in our models we do not make any assumption of alignment of the 
VEVs of the Higgs doublets, viz. we require neither any equality among the Vj nor that 
any of them vanishes. On the other hand, we do require a strong hierarchy of the VEVs; 
indeed, it follows from equation (|16bl) that m e = a/3 \hvi\, etc., and, therefore, 

M : \v 2 \ ■ N = m e : m M : m T . (17) 

2.1 Model I 

In this model we use the type I seesaw mechanism [T2j. The lepton sector contains 
three right-handed neutral gauge singlets Ur and the scalar sector contains a fourth Higgs 
doublet, (j) v . Additionally, there are three complex scalar singlets, which we subsume 
under the symbol r\. The Higgs doublet (fi u is invariant under S, T, and T', while vr and 
7] transform in the same way as D^. A summary of the multiplets and their transformation 
properties is presented in table [H Notice that the symmetries T and T 7 together generate 
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Table 1: Multiplets of model I and their transformation properties. 

a group A(27) under which D^, £r, vr, and r] are identical triplets while (f> decomposes 
into non-equivalent sing letsQ 

1 We may speculate about the full symmetry group of the model. Taking into account also iS, the 
column with caption (.r of table Q] suggests AC to be a generator of the symmetry group; it has the 
sixth root of unity, — ui 2 , in its diagonal. So we might naively guess A(108) = A(3 x 6 2 ) to be the 
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In this way we obtain the neutrino Yukawa couplings 

CP = -y v D L (j) v u R (18a) 
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2 




+y' (v£ R C 1 V3 R T]x + vl R C l v\ R j]- 1 + v[ R C 1 u 2 Rm) + H.c. (18c) 

The neutrino Dirac mass matrix Mp, which originates in (]18aj) . is proportional to the 
unit matrix. We assume the VEVs (i]j)o = Sj to be at a high (seesaw) scale. Therefore, 
the inverse neutrino mass matrix has the form 



AC 1 = I c (b a I , (19) 
\ b a (c J 

with (* = y/y'. The inverse mass matrix ( fl9l) has the typical form of mass matrices 
in renormalizable models based on the group A(27) [14J. This is understandable since, 
as we have pointed out, the fields Dl, u r , and r] behave under T and T' as irreducible 
three-dimensional representations of A(27). 

As explained in the introduction, we need the neutrino mass matrix Ai u to be real. 
The first step in this direction is to assume in the Lagrangian a CP symmetry which 
renders h, y u , y, and y' real. That symmetry is given by 

cp . f Dl iCDi, t R iSCt R , v R zCV* , 
\ -> S0*, 0, 0^, 77 77*, 

where C is the charge-conjugation matrix in Dirac space while 

/I 0\ 

5=001 (21) 




acts in flavour space. This matrix 5 is needed in order to interchange the terms 2 ^2i? 
and 03^3_r in each of equations (EJ). 

Next we discuss the scalar potential of the three rjj. It has six terms compatible with 
T, T', and the CP symmetry: 

3 

V v = E H^'l 2 + Al l^'l 4 ) + A 2 (|^i^2| 2 + |r/ir/ 3 | 2 + |r/ 2 r/ 3 | 2 ) 
j'=i 

+Mi (r/ir/ 2 773 + H.c.) + M 2 (77? + 77I + r/f + H.c.) 

+-^3 (vfv2V3 + vfvm + vfvm + H.c.) . (22) 



symmetry group. However, beyond the symmetries listed in table [T] the model possesses an accidental 
2 <-> 3 interchange symmetry: 

D 2 L -H- D 3 L, ^2fl. O ^3fl, 02 03, ^2fl O ^3i?., ?72 "H- 773. 

Therefore, the full symmetry group is A(216) = A(6 x 6 2 ). For a discussion of the groups A(6n 2 ) see [13]. 
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All six constants in V v are real: /i, Ai, and A2 are real because the potential is Hermitian 
and Mi, M 2 , and A3 are real because of the CP symmetry ( 120]) . The latter three terms in 
the potential are responsible for the relative phases of the VEVs Sj. If we choose negative 
Mi, M2, and A3, then at the minimum of the potential those VEVs will have phases given 
by US] 

args i = cA (23) 

where the pj are integers such that Pi + P2 + P3 = mod 3. This is precisely what is 
needed in order for the matrix Ai^ 1 to be real apart from unphysical phasesH 

It can be shown that the potential (122]) is rich enough to allow for the \sj\ to be all 
different, as needed in our model. 

Since the model has four Higgs doublets, one might be tempted to argue that, after 
switching to a basis in the space of the Higgs doublets where only one of them has a 
non-vanishing VEV, precisely that doublet with VEV corresponds to the Higgs doublet 
of the Standard Model and all other doublets can be made heavy [16J. However, this 
argument is only applicable in the general case where the Higgs potential is not restricted 
by family symmetries. But, if one accepts the possibility of soft breaking of the CP and 
family symmetries through terms of dimension two in the scalar potential, then one can 
apply the above argument in a modified way. Ignoring, for the time being, the gauge 
singlets rjj, we may assume that the term 4>l<j) u has negative sign while the 3x3 matrix of 
mass-squared terms for the <pj is positive definite. Then, the VEVs Vj are induced by the 
VEV of 4> u through the soft-breaking terms [UJ . In this setting, the role of the Higgs 
doublet of the Standard Model is played by <p v , while the <pj are additional doublets that 
can be made heavy. When one includes in the scalar potential terms containing products 
of both singlets rj and doublets 0, the potential becomes even more versatile with respect 
to our goal. 

The models in this paper are designed for the lepton sector. It is straightforward, 
however, to accommodate the quarks by using the doublet 4> v to give them masses, thus 
playing the role of the Standard Model's sole Higgs doublet. In this setting, the doublets 
4>j do not play any role in the quark sector. A problem arises, however, since we have in the 
model a CP symmetry and we know CP to be violated in the hadron sector. This problem 
may be solved by adding to the model one or more extra scalar doublets transforming 
under the various symmetries in exactly the same way as 4> v ; all those doublets will have 
Yukawa couplings to the quarks and, if their VEVs acquire relative phases through the 
mechanism of spontaneous CP violation, then the quark mixing matrix will be complex, 
yet the predictions for the lepton sector will stay unchanged, because the Dirac mass 
matrix Mp will just acquire an overall phase which may be rotated away. 

2.2 Model II 

In this model we use the type II seesaw mechanism [18] . The lepton sector is identical with 
the one of the Standard Model, i.e. no right-handed neutrino singlets are present. The 

2 Clearly, there is a ninefold degeneracy of the minimum (|23[) . hence there will be domain walls. One 
could avoid them through an appropriate soft breaking, for instance a term — Re (771 + 7/2 +773), which 
breaks T but conserves T' and would render the case with real VEVs the deepest minimum. 
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scalar sector contains, besides the three Higgs doublets in 0, three gauge- SU (2) triplets 
with weak hypercharge 1, which we subsume under the symbol A: 



A " V A -A+/V2 J' { } 

We assume A to transform under S, T, and T' in exactly the same way as D^. This 
allows us to write down the Yukawa couplings [19] 



4 A) = f >>.!//• ; A,/;,, (25a) 
+y" (D^C^eAtDsL + Dj L C- 1 sA 2 D 1L + D^C^eA^) + H.c, (25b) 



where 



1 
-1 



(26) 



acts in gauge-SU{2) space. When the neutral components of A acquire VEVs (A°) Q = 5j 
we obtain 

/ yh y'h y'h \ 
M u = y>6 3 y5 2 $8 X . (27) 

V y'fo y'Si y8 3 J 

This is of the same form as the in equation ( fl9lh 

From the transformation properties of the scalar doublets and triplets A, it is obvious 
that the scalar potential cannot have a trilinear term of the form <ft Aecf)* invariant under 
S, T and T'. Absence of such a term leads to a Goldstone boson since the scalar potential 
becomes invariant under separate phase transformations of the <fi and A. In order to avoid 
the Goldstone boson one must resort to soft breaking of the flavour symmetries. We shall 
not pursue this issue further and in the following we shall simply assume that there is a 
satisfying solution which leads, moreover, to real VEVs at the minimum of the potential. 

3 Predictions for neutrino masses and lepton mixing 

We first discuss the predictions of model I. It is convenient to use the weak basis where 
the charged-lepton mass matrix is diagonal. In that basis the neutrino mass matrix is 

M* 1 = VjMJJt. (28) 

Since in our case Ue = U u , we find 

(a c b 

M^' 1 = UlM~ l U: =[ c C~b a |, (29) 

b a (c 
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where 



c-i 

3 

C-i 

3 

C-i 

3 

C + 2 

C-i' 



(a + b + c), 
(a + tu 2 b + cue) 
(a + tub + cu 2 c) 



We know that a, b, c, and ( are real. Therefore, a and ( are real too, and c 
Therefore, the inverse mass matrix (1291) has the form 




with x and w real. 



(30a) 
(30b) 
(30c) 
(30d) 

= b*. 
(31) 



This type of matrices was discussed in [8], where it was shown that its diagonalizing 
unitary matrix, which in the present case is the complex conjugate of the lepton mixing 
matrix, i.e. U*, is of the form [7J [8] 



U 



Ui u 2 u 3 
Wi w 2 w 3 

W{ W* 2 



(32) 



where the Uj are real. Note that the phase convention inherent in equation (1321) is that 



U diagonalizes Ail up to arbitrary signs of the masses. Thus, we have 

WM^U* = diag ( — , — , — J , with ^ = e jmj , ej = ±1. 



(33) 



The phase factors appearing in the effective neutrino mass of neutrinoless double-beta 
decay are identical with the e,-, i.e. 



3=1 



3 

E 



(34) 



There are nine physical quantities in the neutrino sector: three neutrino masses, three 
mixing angles, one Dirac-type CP-violating phase, and two Majorana phases. The mass 
matrix ( |29l) has four real parameters, consequently the model must make five predictions: 
the atmospheric mixing angle 823 is 7r/4, the Dirac-type phase 5 is ±vr/2, each of the 
two Majorana phases is either or it, and the remaining prediction can be gathered from 
inspection of equation (|29|) . whence we deduce that 



MP' 1 ) (m^- 1 ) 



13 



Cab 



M^~ l ) 22 (MM' 1 ) 



23 



(35) 
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This identity is not meaningful in its phase, because M.^ may be rephased at will 



-i 



M M 1 



(36) 



However, the equality of the moduli of both sides of equation ( 135]) is physically meaningful. 
The physical content of this relation will be worked out in the following. 
Due to equations ( |32l) and ( l33l) one has 



-i 



-i 



ii 



13 
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2.3 



3 U 2 

E^> 

i=i ^ 

U j W j 

J= i r 
3 w 2 

E^> 

3 | ,2 

E^- 

i=i ^ 



(37a) 
(37b) 
(37c) 
(37d) 



Notice that (m { ™ ] and (M" are real, c/. equation f l3~T|) . The equality of the 



moduli of both sides of equation fl35l) then reads 



3 «?' 



E? 



E 



3 ^ 



At this point we have to exploit the unitarity condition of U 

UjUji + 2 Re (wjW*,) = Sjji. 
We can transform equation (I38I) into 



(3? 



(39) 



2C/^, 



j<7 
2 



^ 2m 2 

7 J 



E 



x 



(i - >p -i + Uj + ^ + Mi 

1-Uj- U r + / y > 



i<i' 



^ (1 - ^ 



4m 2 



+E 

i<i' 



2/ij/ij' 



(40) 



where we have defined C/j = u 2 . Taking into account that the first line of U is related to 
the mixing angles through 



U 3 

ih 



sin 2 6 



13, 



cos 2 6*i3 sin 2 #12, 

COS 2 6*13 COS 2 #12, 



(41a) 
(41b) 
(41c) 
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the fifth prediction of our model, embodied in equation (140 p . is amenable to numerical 
analysis. 

In the case of model II, the predictions for lepton mixing and for the Majorana phases, 
listed in the paragraph after equation ( )34|) . hold true as well. As for relation ( )40l) . we 
must make the replacement fij — > l/fij in order to obtain the corresponding relation for 
model II. 

In order to evaluate equation (HOI) numerically, we use as input the 2cr ranges of 
Amjj = m 2 — m 2 , Am^ = m\ — m 2 , sin 2 # 12 , and sin 2 # 13 taken from table 1 of |3j. We 
allow the lightest neutrino mass, mojj to lie in between zero and 0.3 eV; this is inspired 
by the extant cosmological bounds on the sum of the light neutrino masses [9]. These five 
parameters, in the specified ranges, form our parameter space. In order to find points in 
this parameter space which are compatible with equation P0|) . we define a figure-of- merit 
function F = \R — L\j \R + L\ + $, where R and L are the expressions in the right-hand 
and left-hand sides, respectively, of equation P0|) . and $ is a function which has the value 
zero if all the parameters lie in the ranges specified above and 10 6 if at least one parameter 
is outside its range. We minimize F and declare a point to be allowed whenever F < 10~ 9 . 
For the minimization we employ the Nelder-Mead algorithm, i.e. the downhill simplex 
method [20] In order to produce the scatter plots in figures [THU for every possible sign 
choice^ of the masses fij — see equation (!33j) — we have run the Nelder-Mead algorithm 
with 10 5 randomly chosen simplices in the parameter space. 

Figures Q] and [2] are for model I. In figure Q] we have plotted the lightest neutrino 
mass mo as a function of sin 2 9 12 ■ Here and in all other figures, the blue (dark) colour 
corresponds to a normal and the red (light) colour to an inverted neutrino mass spectrum. 
The different bands in figured] are associated with different sign choices 6j. These bands 
correspond to the allowed spots and lines in figure |2] which can be deduced from the 
corresponding ranges of the smallest neutrino mass. Let us consider figure [2] for a detailed 
explanation of the sign choices associated with the allowed ranges of m . Of the two 
spots in the left of that figure, the upper one corresponds to (+ + +) and the lower 
one to (+ H — ). In the middle of the scatter plot, both the upper and the lower stroke 
correspond to signs (+ H — ). On the right part of the figure, the upper line corresponds 

to both (H h) and (H ) but the second sign choice holds only in its upper three 

quarters; the lower line, too, is generated by two sign choices: (H ) holds along the 

whole line and (H h) in its upper half. 

Figures [3] and H] refer to model II. The major difference between models I and II is that 
the two lowest bands for mo disappear in model II: there is no band below mo = 10 -2 eV 
in figure [3] and there are no spots in the left side of figure H] This also means that the 
sign choice (+ + +) is not allowed in model II. Otherwise, the interpretation with respect 

3 Note that mo = mi for a normal neutrino mass spectrum and mo = 7713 for an inverted neutrino 
mass spectrum. 

4 This is the method that we have used for producing the scatter plots displayed in this paper. However, 
it is also possible to treat equation (|40l) exactly, since that equation produces, when using as input the 
neutrino masses and C/3, a quartic equation for [7 2 (one must use Ux = 1 — U2 — U3), which is solvable 
through an exact algorithm. We have used this exact method to confirm the numerical results presented 
in this paper. 

5 Note that the sign choices (ei,£27 e 3) and (— e\, — e^, — £3) are equivalent, as can be read off from 
equation ([3D]). 
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■| g-3 I | | | | | 

0.29 0.3 0.31 0.32 0.33 0.34 0.35 

sin 2 8 12 

Figure 1: Lightest neutrino mass as a function of sin 2 9 12 in the case of model I. Here and 
in the following figures, dark (blue) and light (red) colour indicate the allowed range for 
the normal and inverted ordering of the neutrino mass spectrum, respectively. 




,0 



Lightest neutrino mass [eV] 

Figure 2: Effective Majorana mass rripp versus lightest neutrino mass in the case of 
model I. 
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sin 2 6 12 



Figure 3: Lightest neutrino mass as a function of sin 2 12 in the case of model II. 




Lightest neutrino mass [eV] 

Figure 4: Effective Majorana mass mpp versus lightest neutrino mass in the case 
model II. 
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to the signs 6j is the same in both models. 

In both models, the dependence of the allowed range of mo on sin 2 # 13 is very faint and 
does not show up significantly in a plot. For this reason we refrain from showing those 
plots here. 

4 Conclusions 

The recent experimental results on neutrino oscillations have shown that the reactor 
mixing angle #13 is not as small as previously thought. This disagrees with the standard 
version of /i-r interchange symmetry, but not with an alternative version which predicts 
a maximal atmospheric mixing angle and maximal CP violation in neutrino mixing while 
leaving #13 arbitrary. In this paper we have introduced a novel mechanism which realizes 
this scenario. Our mechanism needs a left-handed diagonalization matrix Ui — U u in 
the charged-lepton sector — see equation (j3J) — where U u is the well-known maximal-mixing 
unitary matrix which also appears in ordinary A 4 -based models; our mechanism moreover 
needs a real neutrino mass matrix. We have constructed two models, one based on the 
type I and the other one based on the type II seesaw mechanism. Since with regard to 
the neutrino sector the symmetry structure of our models bears resemblance with models 
based on A (27), we have obtained an additional constraint on the neutrino mass and 
mixing parameters, which can be approximately interpreted as a determination of the 
absolute mass scale of the neutrinos in terms of the mass-squared differences and of the 
solar mixing angle — see figures 1 and 3. Because several sign choices are possible for the 
neutrino masses — see equation (1331) — this determination is, however, ambiguous. Since 
in this alternative version of fi—r interchange symmetry the Majorana phases are either 
zero or ir, there is a rather stringent correlation between the neutrino mass scale and 
mpp, the effective Majorana mass in neutrinoless double-beta decay — see figures 2 and 4. 
It should be emphasized that mpp may in our model assume quite large values, of order 
0.1 eV, which might render it observable in upcoming experiments — see [21] and references 
therein. 

Besides the predictions for the lepton sector, the models have interesting features. 
They do not require any VEV alignment, only reality of the VEVs of the scalar gauge 
singlets in model I / gauge triplets in model II is essential. This reality is easily obtained 
at least in the case of model I, as we have shown here and also earlier in [22]. Possible 
phases of the VEVs of the scalar doublets are irrelevant for CP violation in our mod- 
els. Nevertheless, the models do feature CP violation in lepton mixing, which is indeed 
maximal, i.e. 5 = ±7r/2. The reason is that here the CP symmetry (1201) is not broken 
by complex VEVs but rather by |t>2| 7^ 1^3 1 j which is responsible for the muon and tau 
masses being different, (m^ ^ m T ) [8], [23]. Indeed, in our models the different masses of 
the charged leptons are not obtained from different Yukawa coupling constants but rather 
from different VEVs. This idea has been used before in several other models [211 EE 126] ; 
here we have followed |25J in this respect. 
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